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SUMMARY 

theoretioel  investigation  is  mad©  of  a  circular  jet  issuing 
from  an  orifice  into  an  air  stream  flowing  parallel  to  the  jet  axis 
and  in  the  same  direction  as  the  jet.   The  analysis  is  divided  into  two 
portions  i 

1.  The  region  immediately  downstream  of  the  orifice  in 
which  exists  the  potential  cone  of  uniform  velocity* 

2.  The  region  downstream  of  the  end  of  the  potential  cone 
in  which  the  jet  is  fully  developed. 

The  solution  for  these  two  portions  are  fitted  together  at  the  section 
at  which  the  potential  cone  disappears  jiving  a  complete  solution  for 
the  whole  jet* 

Numerical  computations  are  made  for  several  ratios  of  free  stream 
velooity  to  jet  velocity  to  determine  the  spread  of  the  jet,  velocity  on 
the  axis  of  the  jet,  and  spread  of  the  surface  on  whioh  the  velooity  is 
one  half  the  sum  of  the  jet  velocity  and  free  stream  velocity. 

Squire  and  Trouncer's  results  are  plotted  as  a  comparison  with  those 
obtained  in  this  analysis* 


I.   IHTRODUCTIOB 

The  purpose  of  this  analysis  is  to  provide  a  comparison  of  the 
characteristics  of  a  round  turbulent  jet  discharging  ir.to  a  moving 
air  stream  with  the  known  characteristics  of  a  jet  discharging  into 
air  at  rest,   uoth  jet  and  free  stream  are  assumed  to  heve  the  same 
density  and  the  flow  is  assumed  to  be  incompressible. 

At  the  present  time  there  exists  no  satisfactory  theories  for 
turbulent  mixing.  However,  the  integrated  forms  of  the  equations  of 
motion,  along  with  tome  dimensional  reasoning,  oan  be  used  to  give 
some  useful  results. 

The  solution  of  the  round  jet  in  a  general  stream  has  been 
already  carried  out  by  Squire  and  Trouncer.   In  their  solution  the  shear 

(5) 

stress  is  based  on  Prandtl's  momentum  transfer  theory  which  is  known 
to  be  in  error.   It  therefore  seems  worth  while  to  investigate  this 
problem  using  the  integral  equations  and  also  making  use  of  experimental 
data  for  obtaining  the  shear  stress  in  the  mixing  re -ions. 


Use  of  the  integral  form  of  the  equations  was  first  proposed  by 


(1)  (4) 

von  Rarman  j  one  of  the  earliest  applications  was  by  Polhausen  ,  in 


investigating  laminar  boundary  layers  in  a  pressure  gradient.   Sutton 
originally  used  the  integrated  form  of  the  mechanicsl  energy  equation 

(6) 
in  laminar  boundary  layer  investigations,   uiepmann  and  Laufor  have 

also  applied  this  method  to  some  simple  two  dimensional  turbulent 

flows. 

In  the  integrated  form  of  the  equations  of  motion  the  only 

assumption  nec:3sary  is  that  for  the  velocity  profile  across  the 

mixing  region.   This  oan  be  made  quite  accurately  and  for  e  closer 


comparison  to  the  work  of  Squire  and  Trouncer  the  same  velocity  profile 
that  they  use  will  be  assumed  in  this  solution. 

In  transforming  the  momentum  equation  to  the  integrated  form 
the  term  oonteiainn  the  shear  drops  out.   It,  is  therefore  necessary  to 
introduce  the  mechanical  enerpy  equation.   This  equation,  while  not 
independent  of  the  momentum  equation,  is  necessary  in  order  to  introduoe 
a  term  containing  the  shear  stress. 


The  following  notation  will  oe  used. 


Nozzle  Exit 


End  of 
Potential  Cone 


Fig.  1 
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Distance  from  nozzle  exit  along  jet  axis 

Distanoe  from  jet  axis 

Distance  from  nozzle  exit  to  end  of  the  potential  cone 

Radius  of  the  potential  cone  at  any  point 

Radius  of  the  outer  jet  boundary  At  any  point 

-ias^eter  of  nozzle  exit 

Velooity  parallel  to  jet  axis  at  any  point 

Jet  Telocity  at  nozzle  exit 

Free  stream  velocity 

Velocity  on  the  jet  axis  at  x  — x0 

Velocity  of  the  stream  normal  to  the  jet  axis 

Stream  density 
Shear  stress 


II.   ANALYSIS 

A.   region  between  the  exit  and  the  end  or  the  potential  cone. 
The  three  fundamental  equations  are 


Continuity 


ao-u) 

3X 


,     .dCry) 


ii  omentum 


eu|^,eir|M 


=  o 


I     Q_ 

r  ar 


(rv) 


Mechanical  energy 


eu2^L    +  QWS  |£-  u  |Wrr) 


ax- 


ar 


r    ar 


These  three  equations  are  integrated  between  the 
boundaries  of  the  mixing  region  and  are  then 


(  ontinuity 

4-  urdr  +irbb +^a  ^  -ubb^  =  o 


Momentum 


A. 


u2rdr  +UuiTub 


vol 


ir  i-ub7/bD  +  ilia  ij  -  u;d  ^  =0 


U  db  _ 


liechanical  energy 


1  dx: 

2  JK 


pb 

ur 


dr+K^u*adK-K^-g 


n^ 


rrfNr 


dr 


OL 


A  function  is  assumed  for  the  velocity  profile 
across  the  jet  such  that 


u  =  il  + 


I  -  COS.TT^] 


or 


u  =  Ub  +  u^  |(^ 


where 


b-a 


A 


Substituting;  into  the  continuity  equation 

'  [ub+ti^f^][b-(b-a)d^][4-a)d?] 


ol~>c       Ubk   dx-     "^ 


+Vbb-+Uj3 


Expanding 

dx 


Ub(ab-b2)Jd^  ■ +iUa-b)jV?  +  ^Ydb-b^^)^ 


+^(0-^1^]  +?Tkb  +u.a&  -  ufcb  $  =0 


Integrating; 


d 

35 


-lUab-b1) -  f(a-bf+  **£* (ab-fc?)!, 


+  ^fo-bnt]+TJbb+U«a^  -ubb^  =0 


where 


Differentiating  and  collecting;  terms 


[AQ  +Bb]  4l  +  [Bo  +  Cb]  £2.  +T/bb  =0 


where 


A  =  (Uc-ubXI2+l) 

C  =  (ua-ub)(i2-l,) 

Substituting  into  the  momentum  equation 


LiTub  +Uo!a  4^  -Uub 


•b^b 


d* 


^rr  ~° 


(0 


Expanding  and  integrating 

jLf-uJ&b-b'-;  --^Va-b)  +ubca<t-ukXab-bi;i, 
+  ubvbb  +aia  £f    -  <b  g.  =0. 


■where  n 

differentiating  and  collectinp  terms 


[Da  +  Eb]  ^  +[Ea+Fb]|^-   i-ubVbb  =0 


(2) 


where 

D  =  u* -u|  +  £»  Cufc-ufe)t  +  2IZ  Ubfu.i.-uj 
E.  =  C^1  -  ^)(^.-«J  +CI.-2I  jufcflu-u  J 

The  shearing  stress  now  appears  in  the  mechanical 
energy  equation  and  an  assumption  will  be  made  for  the 
shear  profile.   The  shear  stress  is  a  variable  and  could 
be  determined  from  the  velocity  profile  .  However,  s 
shear  profile  derived  from  an  already  assumed  velocity 
profile  would  have  &  creator  degree  of  error  due  to 
small  errors  in  the  assumed  velocity  profile.   For 
simplicity  an  average  shear  profile  will  be  assumed 
for  the  whole  region. 

From  experiments  with  a  jet  discharging  into  still 
air  shear  profiles  have  been  found  as  shown  in  Fig.  2. 
At  the  jet  exit  the  shear  profile  is  symmetrical.   Past 
the  end  of  the  potential  cone  the  maximum  shear  occurs 
at  about  7?  =.067  to  .75  and  the  shear  profile  has  a  steep 
8 lope  on  the  high  velocity  side.   An  average  shear  profile 
for  the  potential  cone  region  is  shown  in  Fig.  2  and  is 
given  by  


.'-•■  i 


which  is  assumed  to  be  of  the  form 


where 


g^>  = 


u'v' 


In  particular  it  will  be  assumed  here  that  the 
function  g(fe)  can  be  represented  by 


Substituting  into  the  mechanical  energy  equation 


2  3* 


i 


-t^" 


b-4^- 


-_  J. 

~    e 


3*   z^b  "  irv 


Expanding  and  integrating; 


[-Ub+  1  Ufc  (6^-601,  +  I  "btoo  Ubfl, 


+  -f(u*.-ub)3I! 


JL 


ab-b" 


-  i  «b?  -i- 1  <V<^  -^4;  It 


+  U^trbb  4-u^a^    -«tb-§- 


=  tt  ruo_-ab)3(l-,-lg)b  +  it  (u^-ub)3I8a 


where 


I.=    -f>) 


i,  =  %w 

I8=     ^Uyqty 

Differentiating  and  collecting;  term* 


[Ga+Hb] 


da 


+ 


Hcn-Jb 


^4+ub-Vbb  =  Ka+Lb     (5) 


where 

G  = 

J  = 

K  - 
L- 


+  ip  (u«,-ub)2I6 

TT  fU^-Uw)3  I8 

TT    (CU-Ub)3(l7-I8) 


Eliminating  vb  in  equation  (1;  and  equation  (3) 

Ga  +  Hb]  ^+[Ha+Jb]^  -ut [Aa  +Bb]^r 
-uS[bq+Cb]-££-  =  Ka  +  Lb 


10 


b- 


where 


u, 


da 


4- 


+[j-u:c 


H-u 


5][b 


dt  *QW 


^    =KQ+Lb 


[<xa2  +  2@ab  +  ?rb2]    =  Ka  +  Lb 


oc  =  ^A   ?   (3=  Jbi^B.,  T=iui<L 

Using  the  chain  rule   of  differentiating        sLD  -  ^(  )    .  db 

d  *    ~    d  b       dx 
the   above  equation  becomes 


[aota^+2^a.2(3b^^2^b]^  =  Ka+Lb 
(Ka+Lb)  H  =  (2°ca+2(3b)|f  +2pcL  +  2^b 

Eliminating  vfe  in  equations  (1)  and  (2) 


-ub[Ba  +  Cb]^-  =0 


[d-u,*]^ 


-v- 


E-Ubb   a 


it 


dx 


b 


feU^c] 


(4) 


=  0 


Integrating 

Sa2+2e.ab  +■  $  bx  +Z  =o 


fej 


where 


8=^a,     £,£^   ,     <,. 


F-Uv,C 


=   oonatant  of  integration 


Differentiating  with  respeot  to  b 


da. 

db 


ea  +  ct>b 


fe) 


6b  a-  8a 

Substituting  equation  (6)  into  equation  (4) 

[Ka+Lb]^-(2-+2f3b)(-f|±g)-.2(3a  ,2rb 

2  db  Gb+ 6a 


2  db     '  A'      S^  +  eb       Ka+Lb 


■where 


A*  S(eK-SL) 

X    =  KWU-Ks)+KcK'Kft-Lo<.)+L'-0*e-,S(3) 

From  equation  (5) 


6a+eb=  V_z&  i-bl(€"-<t>6)' 


12 


2  db    '  V-?^bXe^4>6)v    Ka+Lb 


+ 


A3b 


Integrating 


*  -v  u     Az'\/-?6+bx(e"L-^g)        v 


bdb 


Ka+Lb 


Ihl»  expression  oen  be   integrated  and  gives 


2-  ^s-e'1 


+Ai£iLr,/'L_»/c\.c  .-) ^  c,    I      b-icTM-C|-fCcv->-0(c,1-c,b>jVc,  .,r 


where 


Cft  = 


LG-KG 


c,= 


Y*  constant  of  integration 
Y  can  be  found  by  the  boundary  condition  that  x  =0 
whore  b  =  d/2.   Obtaining  the  parameters  x  and  b  in 
dimensionlesa  form  the  result  becomes 


£ 


2A,(^-^VE 


0 


±M  -  &)Y*  S  -  €l) 


*S-e2 


+      2AtV~  ^-(*&-6x)' 


<£&  -e1 


+  A_3c 
c, 


lfLc1+>  ,  £ 


S^iif^-i] 


(7) 

The  constant  of   integration    zi   is   obtained  from  the 
boundary  condition  that     a  =  b  =  d/j      at     x  =  0     so 

K 

•    ,Z  -  -^-6u(u«.-U.b) 

The  only  remaining-  constant  to  be  determined  is  -fe 

which  appears  in  the  shear  term.  From  experimental  data 

(6) 
obtained  by  Liepmann  and  Laufer  it  is  found  that  near  the 

jet  exit    SK'fcWx."   °»0072»   From  experimental  data 

(n) 

obtained  by  Corrsin     it  is  found  that  in  the  fully  developed 
region        9^max     =     °»022»      Both  of  these  results  are  for 
the  case   of  a  jet  discharging   into  still  air  or       ~4*L_    =  q 


14 


Differentiating  our  assumed  function  it   is  found  that 
Qtynw occur*  at  the  point  where   £  =  0.692     and      ©j(£)M4^  =•  0.^4  "fe 
The   oonstant    r?  must  be   given  some  average  value  such  that 
^^HAX.for  ^h*  average   shear  profile  will  be    between     0.0072 

and     0.224. 

(8) 
From  experimental  data     it   is  known  that  the  end  of  the 

potential  cone   is  at  a  point  about  9  radii   downstream  of  the 

exit  of  the  jet  for  the   limiting  oase  where     u?  =0.     In  order 

f°r    "  ~37    ~    9     when    u^j*-0     and     a  =0 

3^^  =  0.017157 

and  'k.    =  0.07B88 


15 


6* 


B.   Region  downstream  of  the  end  of  the  potential  cone. 

At  the  point  downstream  of  the  jet  exit  where  x  =xc 
the  potential  cone  disappears  and  the  variable,  a,  no 
longer  appears.   Downstream  of  this  point  the  velooity  on 
the  jet  axis,  uc  ,  will  become  a  variable  decreasing  with  x. 
The  assumption  for  the  velocity  profile  in  the  region  of 
the  potential  cone  was 

a  =  ub  +  ±^  [i  -  cos  ir  (-££)] 

At     x  -  x .     a  ■=  0     and       u  c  =  u  ^ 


U  =  Uu  +- 


"^  [|  -  C45  TV  (  ^j] 


This  reduoes  to 


a  =  Ufc^-    i^^[l+   tOSTT-^] 


or 


U  -  ULb-r 


Ut- 


^.4(/0 


where 


r 


r    -  byx 
dr   =  bdiuL 


id 


Th«  continuity  equation  in  this  region  then  becomes 

Av    ... 


$■    ubbj/xd^  +  Jis^j-bM-f^in. 


4irbb-ukb^-  =o 


This  can  be  reduced  to 


I.V&+  iLb-^« 


1  v  dK 


^°  ^rr  +^  =0 


where 


19  =  ]o-f  T/x)  d^. 
V  =  ut-ub 

Substituting  in  the  momentum  equation 


fe) 


+• 


1  -J 

I  («c-  tO'V j  /JL-f  \yu)  dyuuj  +  Uwvb  b  -  u  *  b  ^  =0 


t    Which  reduces  to 


KubV+I(0V^+[2I9ub+I10V 


b^+2ubirb=o 


(9) 


where 


iio   "J/J-f  (/Jl)^ 


In  the  meehanlcal  energy  equation  the  sheering  stress  Is 
given  by 

The  function  <$(jx)    must  be  euch  that  the  shear  profile  will 
resemble  that  illustrated  in  Fig.  2.  This  will  be 

i  =  constant 
Substituting  into  the  meohanical  energy  equation 

+  uXb-<jb3b|£  =  -bruc-"b)aJoAg(/x)-fV)^ 

which  reduces  to 

+[HXb  *\  iloUbbv  +i  i„bv*]  4t 


whore 


Iit  =j  yu.  g  (/ul)t  (juO  du. 


Eliminating  ▼,  in  equation  (8)  and  equation  (9) 

[4I9ubV\l^l]4*  +  kl9ub  +  I,.V]b$ 


+  2ub[-I9V&-±I.b^]   -0 


db  .     l.ub+I.V     _dV 

b  2laUv,  +  I,oV  V 


Integrating 


log  b=--g-log[ii9uv,-hllDv]--^loqV-t-log  C, 

b  - ,     c>  oo 

db  _      b      IsUfc-kUV      dV 
6t  ~      V*(ei9ub+I10V)*  dy. 


Eliminating  vfe  in  aquation*  (8)  and  (10)  and  substituting 
for  -jr   above  m    ret   a  single  equation  which  is 

I9I,0Ub  4-I9I„ubV+iI,0I„V;L   dV  _  _  21^  w  % 
C2I9U^I„V)1A  d*  C, 

This  may  be   integrated  to  give 


8I„I 


17.    -»-0 


C,I„  I 


feubX  =  - 


1 


10 


31,. 


V 


Ll 


Vt 


-2 


i, J 


V 


Vx 


'I. 


(±* 


•i0 


V 


2f  Ut4V 


'A 

+-C. 

(,z) 


The  oonstant  of  integration  C  ^  may  be  found  by  the 
boundary  oondition  that  at  x=xc,     Hs.~ub' 


-■Li„    UJJfy,    -  c  x„  j*Uk54t 


where 


t 


U<L-^-^.)Uk 


The  oonstant  of  integration  C,  may  be  evaluated  from 
these  same  boundary  conditions  and  from  equation  (5)  -when 
a  =  0 


Substituting  these  values  of  C,  and  C,  the  final 
solution  is  obtained. 

"^"    l[i,-i,.]^+i(rn-UO-t.) 


% 


A*&:-%^+Wm^. 


X\Zl<u~<fc 


■i- 


k^-% 


(13) 


where 


I 


a*».ub+v 


a, 


V 


Equation  (11)  becomes 


-L^-^][i,.^+(ai.-iJI 


^Mbvuvf^) 


^U  +T   _J  ^I±. 


*1 


xc  may  be  found  by  solving  for  b  when  a=0  in 
equation  (5).   Using  this  value  of  b  find  the  corresponding 
value  of  x  in  equation  (7)  which  will  be  xc  . 

Equation  (13)  gives  x  as  a  function  of   (uc-uj. 
Equation  (14)  jrives  b  as  a  function  of   (u  -uj,   It  does 


not  seem  necessary  to  ^et  an  explicit  expression  for  x 
and  b.   By  assuming  values  of   (ut-u^)  oorrespondinr  values 
of  x  and  b  are  found  and  b  is  plotter,  against  x  in 
Fig.  5. 

Xo>;  the  oonstant  Z  appearing  in  the  shear  stress 
expression  must  be  evaluated.   In  the  fully  developed  region 
downstream  of  the  end  of  the  potential  cone   Q  (/w)*ux-  3.022. 
Prom  the  assumed  function  of  °( C/J-)  it  is  found  that 
9M*ax  =  0.4382   or  I   *  0.05116. 


22 


III.  DISCUSSION  • 

The  results  plotted  for  the  region  containing  the  potential 
oone  are  the  results  obtained  by  Squire  and  Trouncer,   In  using  the 

derived  equation  for  x  and  b  (equation  7)  in  the  potential  oone 

* 
region  the  final  results  obtained  did  not  appeer  physically  reason* 

able,  particularly  at  hidier  velocity  ratios  of  free  stream  to  jet. 

Due  to  the  enormous  amount  of  computation  necessary  to  recheck  these 

results  Squire  and  Trouncer's  results  were  used  in  this  region. 

The  very  close  comparison  of  the  two  analyses  in  the  fully 
developed  region  would  seem  to  indioate  that  the  integrated  equation 
solution  that  has  been  used  here  should  give  fair  agreement  in  the 
potential  cone  region. 

The  basic  analysis  of  the  potential  oone  region  is  known  to  be 
correct  through  equation  (5),   This  equation  gives  a  relation  between 
the  spread  of  the  jet  and  the  width  of  the  potential  cone.   The 
remainder  of  the  analysis  within  this  same  region  or  the  computations 
may  be  in  error.  However,  the  metnod  of  solution  is  outlined  for  any 
future  investigation  of  this  problem. 
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AFPEHDIX  I 

Evaluation  of  dimenaionless  integrals  appearing  in  the 
equation*  of  motion. 

The  dimensionlecs  integrals  appearing  in  the  equations 
of  motion  all  may  be  expanded  into  integrals  of  forms  found 
in  standard  integral  tables.  Following  are  the  results  of 
these  integrations. 

I,  =  j,  Updy  =J[i-cosir)Z]ol^  =  -1.00 
I,=  J  yft(y)d>?  --  r[i-co57r»J»d»  =  -o.702 


r° 


I.rjftydyi      =j[)-C05TT^]2d^     =     -1.5-00 
l^l^fyfy  =J%0-COS  TT^]\d>^=     -  IJ55 


Is=p3fy)d^  -J°[i-cos^^]3a^  -  -z.soo 


I9=JoyLL-f(|x)cV   =-Jp+C0^^/\l/^4*A    *    °'Z9Q 


1 ,0 - Jo ^-f  C/a,) dyu.    -  r  [itcosrr^]V^  =  0.346 


1  = 


/J   (ia)d^  -f  Q-+-C0S  TTyU}  ya  dyu.  -    C.4-85 

J  o 


lit.*]  Mgc/jg-f'^d^L 


ftLSintlyujjt  COSTTU.    ^-f-L,  (CCSTryU.-hTTJLA.Sinrr^ 


-0.03982 


APPENDIX  II 

Integrations  necessary  for  solution  of  the  equations. 

A.  Integral  appearing  in  the  equation  for  the  region 
containing  the  potential  cone. 

T-f      bcfb 
I  Ka-t-Lb 

Substituting:  for  a  froc  equation  (5) 


1  = 


bdb 


O 


which  reduces  to 


r 


I  =  c0 


bdb 


bWc.'-cb*' 


where 


c.= 


s 


o 


LS  -K€ 


1         Cls-k<^ 

Rationalising  the  denominator 

f  (bx-  bV^-c^jdb 


I  =c0 


b"(ci+.)-C- 


I.-c. 


bxdb 


b^+ij-c 


~    Co 


1  = 


b^db 


bxcTnFc; 


i- 


f 


b"db 


[b^c.  +  r-c^bVc^i     -he,] 


c, 


+ 


2(c,+  i) 


Z(C^+i) 


[^1  bVo+T-c,        bVc^i'+c, 


db 


£  + 


c, 


b^TcTn-c,       bVci+T+c,  Je(c^+»; 


db 


z^L2b+^^(b^~^  -^log(fcW5?T  <-C|) 


2(ci+.){/D  +  V^T    ,09   bVc^T  +"c: 


-// 


i  = 


Le-t        c,1  -c,b1=  ClR*0 

-2c^bcfb=    ZqRod^o 
bdb=  --R0dR0 


I 


// 


i*        C^-C^p7" 


-I 


(Cx*!)*.1     -'-* 


J 


2  + 


C, 


c, 


R.VcTTT-  ^       ^o-vT^h 


-»- 


c, 


^  J 


dl?, 


B.  Integral  appearing  In  the  equation  for  the  rerlon 
doenstream  of  the  end  of  the  potential  cone. 


I 


f  4f:u;+4|*u,.v  +  vl 


V%(z^ub4-V) 


M. 


dV 


Let 


t    = 


V  = 


dV  = 


V 


*3 

ub+V 

z^ 

:0ubt^ 

1  - 

--t*- 

4- 

r«           _l 

0  -t*) 


"M  "2- 


dt 


iKevn 


1=7 


3T 


A 


£^^-^]tv[i^(^T]t 


dt 


± 


L  = 


( 


HK 


3  i,  t: 


*&Kfe' 


-fc 


I., 


ft 

-j 


h 
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